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I. INTRODUCTION 



Accompanied by the development of microtechnology, such as MEMS and /iTAS, there 
is increasing interest on the methodology to realize a desired motion of a micro object in a 
solution environment. It is well known that the principle to create an electric motor in a 
macro system is not applicable to micro system because of the enhanced sticky interaction 
and higher viscosity in micrometer sized system. On the other hand, living organisms 
generate various motions on microscopic scale under isothermal condition. Despite the 
past intensive studies 1,2 , the underlying mechanism of the biological molecular motors has 
not been fully unveiled yet. Under such development status of science and technology on 
micro- motor sat the present, we report a simple motoring system which work smoothly in a 
microscpic scale. 

Recently, we found that rhythmic motion is generated for an aqueous droplet in an 
oil phase under DC voltage on the order of 50 - 100 V. We have already reported some 
experiments and models for a w/o droplet under DC electric field. 4,1 There are some reports 
about experiments of w/o droplets moving under electrical field 5,6 , bouncing and being 
absorbed on a surface between water and oil , deforming and spliting 8 ' 9 . Manipulating this 
kind of droplet, which is interesting as the model of the cell 10-12 , the micro-sized reactor, 
by optical tweezers 13 , by micro channel 11,10 . And manipulating the cells or micro objects by 
electrical field has been attempted 16 , to know manipulating this kind of droplets in detail 
will help this in the future. In the present article, we will show that rhythmic motion on 
micro-droplet is induced under the DC potential on the order of several volts. We will 
also propose a simple mathematical model to reproduce the rhythmic motion and mode 
bifurcation. 



II. EXPERIMENTAL 

A schematic illustration of the experimental setup is given in FIG.l. A water droplet 
was suspended in mineral oil on a glass slide, and constant voltage was applied to the 
droplet using cone-shaped tungsten electrodes. Droplet motion was observed using an optical 
microscope (KEYENCE , Japan). 

The w/o droplet was generated using a vortex mixer as follows. We prepared mineral 
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oil including surfactant: lO/xm surfactant, dioleylphosphatidylcholine (DOPC) (Japan), was 
solved in mineral oil (Nacalai Tesque, Japan) by 90 min sonication at 50 °C. 2/j.l ultrapure 
water (Millipore, Japan) was added to 100 /d of the prepared mineral oil, and then agitated 
by a vortex mixer for approximately 3 s. 



FIG. 1. Schematic representation of the experimental setup. Mineral oil containing water droplets 
was placed on a glass slide and couple of electrodes was situated inside the oil phase. V: Applied 
DC voltage. L: distance between the electrodes. 

III. RESULTS 

FIG. 2 exemplifies the motion of a droplet under DC electric field, indecating the occur- 
rence of the periodic go-back motion between the electrodes accompanied by the increase 
of the electical potential. In the experiments, we observed two following types of behavior: 
oscillatory and stationary. These behaviors switch each other depending on the applied volt- 
age. When the distance between two electrodes was 213/im [FIG. 2 (a)], the droplet started 
the motion with the applied voltage above 16.3 V. When the distance between two electrodes 
was 141/im [FIG. 2(b)], the droplet started moving with the applied voltage above 13.7 V. 

FIG. 3 shows the diagram of the mode of droplet behavior depending on the applied 
voltage with the size of the droplet is — -. When the distance of two electrodes is below 
approximately 70/im , droplets are sticked to an electrode (adhered). The diagram indecates 
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that the threshold of applied voltage is roughly propotional to the distance between two 
electrodes. 

(a)L = 213nm (b)L=14^m 
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FIG. 2. Spatio-temporal diagram on the motion of a droplet with the diameter of 34 m at (a) 
L=213m, (b) L= 141m. ©Bifurcation from the stationary state into an oscillatory state is induced 
by the increase of the applied voltage. 

IV. DISCUSSION 

We propose a model to describe the oscillatory-stationary motion of w/o droplets. In 
an eqation of motion at a micrometer scale, a viscosity term is more dominant than inertia 
term because the Reynolds number, R e , is rather small; R e = pvd/t] ~ 10~ 9 <C 1, where 
p(~ 10 3 kg/m 3 ) and 10 3 Pas) are the density and viscosity of the mineral oil, respectively, 
and t>(~ 10 _4 m/s) and <i(~ 10~ 5 m) are the velocity and diameter of the water droplet. 
Therefore, an over-damped eqation of motion under the constant electric field, E, is given 
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FIG. 3. Phase diagram for the mode bifurcation between rhythmic motion and stationary state as 
obderved for the droplets with different diameter, where each point represent the threshold value 
on the bifurcation. The two arrows are correspond to the mode bifurcation as shown in FIG. 2 

by 

kx = qE + aVE 2 (1) 

where k(= Qirr/d ~ 10 _7 kg/s) is a coefficient of viscosity resistance, and kx represents the 
viscosity resistance for a moving droplet with diameter d and velocity x. qE and aVE 2 
indicate an electric force and a dielectric force acting on the droplet with charge q and 
polarizability a (II) 1 '. 

Here we assume that the time-dependent rate od the charge, q, is described as 

q = _/3 ex 3 _ 1 (2) 

where /3 is the proportionality coefficient, and e is the constrant is proportional to the 
magnitude of the electrical field. The first term of this means the time-dependence of charge 
is in proportion to the number of lines of electric force. The second term means the charge 
leak, and to is the relaxing time. 

We would like to consider the condition where the droplet stays on the same position 
between two electrodes. When the electric field can be written as E = (E x ,E y ), The force 
on the second term in the right hand of Eq. (1) is caused by the number of the lines of 
electric force penetrating the droplet. Comparing with the size of droplet, the change of E x 
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along x axis can be neglected. By considering the symmetry of the system, we simply adapt 
that the change of E 2 along x axis is written as 

£ 2 = e(-(* + l) 2 (*-l) 2 + 2). (3) 

Then the x component of eq.(l) is given as 

kx = qE x + ad x E 2 
qE x + ad x E 2 y 
=qE x — Aaex(x + l)(x — 1) 
x=fq-^fx(x + l)(x-l) (4) 

For simplicity, we introduce the following parameters: E x /k = a, Aae/k = e and kefi/Aa = 
7. Then Eq.(2) and Eq.(4) can be written as 

x = — ex(x + l)(x — 1) + aq (5) 
q= -jex^-q/to (6) 

FIG. 4 shows the result of numeric calculation with these equations, where the time and 
space scales ,T and X, are arbitrary. The change of the distance between the electorodes 
corresponds to the change of the magnitude of the electric field. For example, if L becomes 
larger, a and e become smaller. In FIG. 4, a and e in (b) are larger than those in (a). The 
frequency of the back-and-force motion of the droplet is faster when the electric field between 
the electrodes becomes stronger. Thus, our numerical model reproduces the essential aspect 
of the rhythmic motion of a droplet under DC voltage. 
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